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We perform detailed analytical and numerical studies of a recently proposed method for a spectroscopic 
measurement of the superfluid fraction of an ultracold atomic gas [N. R. Cooper and Z. Hadzibabic, Phys. 
Rev. Lett. 104, 030401 (2010)]. Previous theoretical work is extended by explicitly including the effects of 
non-zero temperature and interactions, and assessing the quantitative accuracy of the proposed measurement 
for a one-component Bose gas. We show that for suitably chosen experimental parameters the method yields 
an experimentally detectable signal and a sufficiently accurate measurement. This is illustrated by explicitly 
considering two key examples: First, for a weakly interacting three-dimensional Bose gas it reproduces the 
expected result that below the critical temperature the superfluid fraction closely follows the condensate fraction. 
Second, it allows a clear quantitative differentiation of the superfluid and the condensate density in a strongly 
interacting Bose gas. 
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I. INTRODUCTION 

Ultracold atomic gases provide a highly controllable exper- 
imental setting for studies of many-body quantum phenomena 
such as Bose-Einstein condensation [1] and superfluidity J2]- 
4 1 (for a review see Q). The physical phenomena studied 
in these systems are often analogous to those occurring in 
other many-body systems, in particular solid state materials 
and liquid helium. Moreover, the flexibility in experimentally 
designing the geometry and interactions in atomic gases of- 
fers the possibility to experimentally access physical regimes 
which are of theoretical interest but could so far not be ob- 
served. At the same time, the experimental probes used in 
atomic physics are often quite different from those tradition- 
ally used in condensed matter experiments and discussed in 
the theoretical literature. In particular, for atomic gases there 
is no general experimentally established method allowing a 
quantitative measurement of the superfluid fraction, which is 
traditionally defined through the fluid's response to rotation 
and the emergence of a non-classical moment of inertia O. In 
classic experiments on liquid helium [7 1, under rotation of the 
walls of the container a perfect superfluid remains metastable 
in the zero angular momentum state as long as the rotation 
rate does not exceed the critical velocity for a superfluid flow. 
More generally, the fraction of the fluid which does not rotate 
with the walls quantitatively defines the superfluid fraction. 

Only very recently some ideas on how to measure the su- 
perfluid fraction of an ultracold atomic gas have been formu- 
lated [8-10 1 . Specifically, in Ref.[l0]a spectroscopic method 
was proposed, which closely follows the traditional definition 
of the superfluid density but allows a signal readout which 
plays to the strengths of atomic physics. This proposal builds 
on the recent developments in the use of optical fields to in- 
duce artificial gauge fields for ultracold atoms ifTTl . The key 
idea is that if slow rotation of the gas is induced by an optical 
gauge field, this creates a natural coupling between the exter- 
nal (angular momentum) and the internal (hyperfine) atomic 
degrees of freedom. Measuring the populations of hyperfine 
states in an atomic cloud then allows a direct readout of the 
angular momentum induced by the rotation, and thus a mea- 



surement of the moment of inertia and the superfluid fraction. 

The basic connection between the hyperfine populations 
and the superfluid fraction was pointed out in Ref. 10 by con- 
sidering the difference between a perfect superfluid with no 
angular momentum and a fully relaxed gas which exhibits the 
classical value of the moment of inertia. Here we extend this 
theoretical work in several ways. First, we include in our cal- 
culations the effects of non-zero temperature and interactions 
in the gas, which broaden the distribution of angular momenta 
around zero for a metastable superfluid, and around the clas- 
sical value set by the imposed rotation for a fully relaxed gas. 
This allows us to assess the quantitative accuracy of the pro- 
posed measurement, and to estimate the experimental param- 
eters which in practice would allow a good compromise be- 
tween the theoretical accuracy of the method and the experi- 
mentally relevant size of the readout signal. Second, we ex- 
plicitly calculate the expected experimental signal in two im- 
portant cases: For a weakly interacting Bose gas we show that 
the spectroscopically deduced superfluid fraction closely fol- 
lows the condensate fraction below the critical temperature; 
this confirms that the proposed method gives the expected re- 
sult in this well understood limit. On the other hand, in the 
limit of strong interactions the superfluid and the condensate 
fraction of the gas can be quite different, as is known from 
the case of liquid helium [12]. We show that in this limit the 
spectroscopic measurement is sufficiently accurate to allow a 
clear experimental distinction between the two quantities. 

The paper is organized as follows. In Section[ll]we lay out 
the theoretical background on the concept of superfluid den- 
sity and its connection to hyperfine populations in an atomic 
cloud rotated with use of optically induced gauge potentials. 



In Section III we give some more details on a specific im- 
plementation of gauge fields in an atomic system. In Sec- 
tion IV we discuss the quantitative theoretical corrections to 



the mapping between the hyperfine populations and the su- 
perfluid fraction. To illustrate the effect of these corrections 
we first consider a normal gas at non-zero temperature; this 
already allows us to anticipate suitable experimental param- 
eters which lead to sufficiently small theoretical inaccuracies 
and sufficiently large experimental signals. In Section [V] we 



2 



extend our calculations to interacting gases, at both zero and 
non-zero temperature. Our results for the expected experi- 
mental signals in both weakly and strongly interacting gases 
are presented in Section VI Finally, Section VII contains a 
summary of the paper. 



II. SUPERFLUID DENSITY 

The concept of a superfluid density, or superfluid fraction, 
originates in the two-fluid model for the hydrodynamics of 
superfluid 4 He, proposed by Tisza fT3lfT4l and Landau lfT31 . 
The fluid, of total density p, is assumed to consist of a super- 
fluid component, of density p s , which has vanishing viscosity 
and flows without dissipation, and a normal component, of 
density p n = p — p s . Landau proposed |15| how to measure 
these separate components. He envisaged taking superfluid 
helium at rest in its container, and slowly rotating the walls at 
a constant angular velocity uj. The normal component comes 
to rest with the rotating walls; however, the superfluid compo- 
nent is unaffected and remains at rest. Since only the normal 
component moves, the moment of inertia of the fluid is deter- 
mined by p n , and its ratio to the expected classical moment 
of inertia (defined by the total density p) provides the normal 
fraction p n / p and hence the superfluid fraction 1 — p n /p. Note 
that it is necessary that the trap is not perfectly rotationally 
symmetric (i.e. the walls must be rough), so that the normal 
fluid can relax into the steadily rotating state and come into 
equilibrium by changing its angular momentum. 

This method was implemented for superfluid helium in the 
classic experiments by Andronikashvili Q. In those exper- 
iments it was not the container that was rotated, but a stack 
of disks embedded in the fluid. Still, the disks drag just the 
normal fluid, so measurements of the moment of inertia of the 
disks (using a torsional oscillator) allowed a determination of 
the normal and superfluid fractions. 

The non-classical moment of inertia arising from the super- 
fluid component provides the standard definition of the super- 
fluid fraction |6|. To discuss this theoretically, it is customary 
to consider the fluid to be contained in a ring-shaped toroidal 
vessel with a radius R much larger than its transverse dimen- 
sions AR, cf. Fig. [T] In this case, the classical moment of 
inertia for N atoms of mass M is given by I c \ — NMR 2 . We 
shall assume this geometry throughout this paper — in part for 
theoretical simplicity, but also for practical reasons discussed 
further below. The superfluid fraction is then defined 1 6 1 by 
the average angular momentum (L) picked up under rotation 
of the vessel with an angular frequency u: 
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The limit ui — ► of slow rotation of the vessel is required 
such that the velocity of the walls of the container, ujR, does 
not exceed the critical velocity of the superfluid, ujR < u C rit- 
Since finite-size effects can be much more important in ultra- 
cold atomic gases than in typical experiments on superfluid 
helium, it is helpful to expand further on this point. For a 
very small rotation rate, namely when ui < H/2MR 2 , even 




FIG. 1: Geometry of the vessel considered in this paper: a torus of 
radius R with transverse dimensions AR <C R. is the gauge 
potential introduced by light fields, corresponding to rotation with 
angular frequency oj c ff. This induces a superfluid flow with speed v s 
in the direction opposite to A v . 



an ideal Bose gas shows a non-classical moment of inertia, 
and could thus be considered to be superfluid [17|. (For 
uj < h/2M R 2 , the lowest-energy single particle state remains 
the state with vanishing angular momentum, so an ideal con- 
densate has L = 0.) We define superfluidity in the strongest 
sense of Ref.[l2l which is the sense that is conventionally ap- 
plied: the angular momentum is measured with an imposed 
rotation frequency in the range U/2MR 2 < u < v a i t /R. The 
lower limit excludes the ideal Bose gas as a superfluid. For a 
large system, R ^> h/Mv cl -i t , this constitutes a wide range of 
frequencies. For a weakly interacting Bose gas, v cr a ~ h/M^ 
(with £ the healing length), and so the range of frequencies is 
wide for R > £. 

This method could, in principle, be applied to ultracold 
atomic gases, using a rotating deformation of the trap to re- 
present the rotating walls of the container. However, in prac- 
tice this is difficult, as a measurement of the induced angu- 
lar momentum or mass flow is difficult for ultracold gases. 
For harmonically trapped gases, an ingenious method of mea- 
suring the angular momentum has been applied |[8l [181 \19\ . 
A theoretical proposal has shown how the superfluid density 
could be extracted more generally if local imaging is possi- 
ble 1 9 1. In this paper we explore an alternative theoretical pro- 
posal [10 1 in which the rotation is simulated by an optically 
induced gauge field. A key feature of this method is that it 
allows one to measure the angular momentum spectroscopi- 
cally and hence deduce the superfluid fraction of an ultracold 
atomic gas. 

The basis of the idea is the coupling of light with orbital an- 
gular momentum lfl6l to internal atomic spin states, thereby 
creating an azimuthal gauge field A v [20|. The azimuthal 
gauge leads to the same effects as rotation, albeit in a slightly 
different manner. In the presence of the gauge field one 
must make a distinction between the "canonical" momentum 
Pcan = — iTiV and the "kinetic" momentum p = p can — A. 
In the absence of a gauge field, a superfluid that is at rest in 
the toroidal vessel has no winding number of its phase (no 
vortices) and corresponds to the case of vanishing canonical 
momentum p can = 0. This does not change when the gauge 
field is introduced. However, as the gauge field is turned on, 
the superfluid picks up a non-zero velocity p/M = —A/M. 
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On the other hand, for a normal fluid, as the gauge field is 
switched on, the fluid will always stay at rest with the walls 
of the container (provided they are rough). Thus, compared 
to the rotating container discussed above, here it is the super- 
fluid that moves while the normal fluid stays at rest. The case 
of the gauge field is, in fact, exactly equivalent to a rotating 
vessel, but where one views the system in the rotating frame 
of reference and so experiences the trap (and normal fluid) to 
be at rest. 

To make this discussion more precise, note that when the 
optically induced gauge field is on, the atoms experience an 
effective dispersion relation which can be written as 



one can deduce the angular momentum expectation value 



E ~ En 
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where I is the angular momentum in units of h, such that it is 
quantized to integer values. M* is a new effective mass of the 
atoms, and £* is the rotational shift due to the gauge field. (We 
will derive this in the next section.) This effective dispersion 
is equivalent to an azimuthal gauge field A v — M* /R. It 
corresponds to being in a frame of reference rotating with an 
effective angular velocity 
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e.g. a particle with £ = will have an angular group velocity 
(l/h)dE/d£ = — w e ff . Following the above discussion, if w e g 
is slowly increased from zero, the superfluid will remain in 
its original state with (L) = but will flow with speed v s — 
ui e ffR in the direction opposite to A v , whereas the normal fluid 
will pick up an angular momentum of hi* per particle but will 
stay at rest. 

A key feature of the proposal [ 1 1 is that since the gauge 
field is generated by mixing internal hyperfine states of the 
atoms, this provides a natural coupling between internal and 
external degrees of freedom. Spectroscopically measuring the 
population of the different internal spin states allows the aver- 
age angular momentum per atom h(£) to be deduced. We will 
focus on the population difference of the hyperfine states |+1) 
and | — 1) in a three-level system, but this could be adapted to 
other internal structures. (The three-level system will be de- 
scribed in further detail in the next section.) For a single atom 
with angular momentum £ we define the population difference 
as 



A^^iv-iwr-iv+iwr 



(4) 



A measurement of the populations A^ + i and N-\ for a gas of 
such atoms leads to the fractional population difference Ap, 
which may be expressed in terms of Api as 
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Within the assumption that 

Apt ~ Ap 



can be expanded to first order, 

A P '£ + 0{£ 2 ), (6) 



Nh " W " Y,Mi) ~ Ap' ■ 
Putting this back into Eq. ([TJ, one gets 
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(8) 



where the appropriate moment of inertia I c \ = NM*R 2 has 
been used. For a perfect superfluid we would expect to find 
Ap = Apo, whereas for a normal fluid we would expect to 
find Ap = Apo + Ap'£* , thus allowing us to distinguish be- 
tween the two. 

The main goal of this paper will be to analyse the quanti- 
tative accuracy of ([8]), allowing for corrections that arise from 
the higher-order terms that are neglected in Eqs. (2]6i. How- 
ever, for now, note that the spectroscopic technique should 
show a clear qualitative signature of superfluidity. For a nor- 
mal fluid (or the normal fraction), it does not matter in which 
order one increases I* and cools the gas to its final temper- 
ature: that is, these two operations "commute". However, 
for the superfluid fraction these operations do not commute: 
if one first cools at £* =0 and then imposes non-zero £*, 
the superfluid is pushed into a (metastable) state in which 
it is moving with respect to the walls of the container; if 
one first imposes non-zero £* and then cools, the superfluid 
will be formed at rest with respect to the walls. Thus, de- 
pending on the order, the system is led either to the non- 
relaxed metastable superfluid condensed in a state of vanish- 
ing (canonical) angular momentum £ c = 0, which we label 
"SF", or to the relaxed superfluid, condensed in the ground 
state £ c =£*, which we label "RSF". The two cases will have 
different fractional hyperfine populations, Ap RSF ^ Ap SF , 
so the change in population allows a clear qualitative signa- 
ture of metastable superfluid flow. 



III. OPTICALLY DRESSED STATES 

There are several well-established theoretical proposals for 
how optical fields can be used to create fictitious gauge fields 
in neutral atomic gases ifTTI . The scheme we follow here is 
closely related to that implemented in the experiments of the 
NIST group l2Tll22l . However, it is adapted to generate an az- 
imuthal vector potential, by using optical beams with orbital 
angular momentum |20l . As described above, throughout this 
work we assume that the gas is confined in a toroidal trap 
(Fig.[T]i with radius R large compared to its transverse dimen- 
sions AR. This simplifies the experimental implementation 
of the azimuthal gauge field, as it will be sufficient to require 
that the optical fields are uniform only over the range AR. 

We consider atoms with three hyperfine levels 11211 in their 
electronic ground state, e.g. 23 Na with F = 1. The de- 
generacy of the three hyperfine states Mp = 0, ±1 is lifted 
by applying a weak external magnetic field B, thereby in- 
ducing a Zeeman shift AE = Z ■ Alp with an energy gap 
Z = gpfiBB between the hyperfine states. These states are 
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FIG. 3: Coupled (solid lines) and uncoupled (dashed lines) energy 
FIG. 2: Sketch of a two-photon Raman transition between hypernne bands of a three-level system, from Ref. [lO] Parameters: S = 
levels, adapted from Spielman et al. (2l]|22). 0.5h(A£) 2 /M R 2 . The smooth curves interpolate between the al- 

lowed integer values of I. 



coupled by two co-propagating Laguerre-Gauss beams with 
frequencies u\, L02 and orbital angular momenta t\, £2- The 
frequencies are chosen such that they are detuned from any 
actual electronic transition, thus single-photon transitions are 
suppressed. Instead, the hyperfine states are coupled by two- 
photon Raman transitions, cf. Fig. [2] For every two-photon 
transition the atom experiences a net change in its centre of 
mass angular momentum of A£ = — while the change in 
linear momentum can be neglected. Hence there exists a cou- 
pling between \M F = -1,£ - M) and \M F = 0,1) as well 
as between \M F = Q,£) and \M F = +1,1+ At). The two 
light beams are slightly detuned from the Raman two-photon 
resonance, with detuning S = (cji — CJ2) — Z/Ti. Including the 
kinetic energies of the different angular momentum states and 
applying the rotating wave approximation [21], one arrives at 
the full Hamiltonian, H{£)/h, 



^{i + AlY-5 n R /2 

n R /2 





n R /2 

h {£-A£) 2 + S 



2MB 2 



(9) 

The £ dependence has been made explicit; each atomic state 
is defined by the amplitudes of the three hyperfine states and 
its angular momentum £. Here fl R is the two-photon Rabi 
frequency which describes the coupling between hyperfine 
states. The effect of the Zeeman splitting is given by 8 11231 . 
For each £ there are three energy eigenvalues of ([9]), corre- 
sponding to three dressed energy bands. 

In Fig. [3] we show the results of a numerical diagonaliza- 
tion of the Hamiltonian d9). The three uncoupled energy levels 
for fl R = are shown as dotted lines. When the light is on, 
fl R 7^ 0, these levels are mixed and lead to the energy lev- 
els of the dressed states (solid lines). For non-zero detuning 
S the minimum of the lowest band is displaced to a non-zero 
angular momentum £* . Provided that all atoms are restricted 
to states in the lowest-energy dressed band, the atoms experi- 
ence an effective dispersion relation of the form (j2j) in which 
the non-zero £* plays the role of an azimuthal gauge field. 

Throughout this paper, we assume that only this lowest 
dressed band is occupied. This is justified if the chemical po- 



tential /1 and temperature kT are small compared to the band 
splitting, which is of order Ml R . Hence we require that Q, R 
be sufficiently large. The lowest band will be referred to as 
ejj {£) in SectioniVl To obtain the results shown in Section VI 
we determine ejj {£) by numerical solution of |9]). However, to 
allow an understanding of the general trends, we derive some 
analytic expressions which are valid for large £l R , where the 
bands are far apart from each other and the lowest band is 
nearly parabolic. Perturbation theory in 1 /Vl R shows iflOl that 
the minimum of the dispersion relation is shifted from £ = 
to 



t ~ - V2—A£ + 0(l/n 2 R ) , (10) 
and the bare mass is increased to an effective mass M* , with 



M* ~ M 1 + 



V2h(A£) 2 

MR 2 n R 



■o(i/n*) 



(ii) 



Similarly, perturbative calculations of and 
|^ +1 (^)| 2 show that they have equal and opposite contribu- 
tions linear in £. The difference |-0-i(-Q| 2 — | , 0+i(^)| 2 can 
indeed be written as a series in £, as in Eq. |6|, with 
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0(i/n R ). 



(12) 



(13) 



MR 2 n R 

As parameters representative of experiments, throughout this 
paper we consider sodium with M = 23 m p (where m p is the 
proton mass) in a trap of radius R = 10 /im. A typical value 
of the effective mass is M* rj 1.15M (at a two-photon Rabi 
frequency Q R = 2ir x 100 kHz and for A£ = 50). 



IV. CORRECTIONS 

In the preceding sections, we have summarised the theoret- 
ical proposal of Ref. 10 This showed how to relate the su- 
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perfiuid fraction to a spectroscopically determined hyperfine 
population imbalance [see Eq. dHb]. The quantitative accuracy 
of this relation relies on the validity of the termination of the 
Taylor expansions in Eqs. <|2j and |6]l at quadratic and linear 
orders respectively. If these (terminated) expansions were ex- 
act, the superfluid fraction would be perfectly determined by 
Eq. (IS}. In practice, higher-order corrections do exist, i.e. 



E — -Eparabolic + c£ + . . . 
A Pi = (Ap^unear + c' £ 2 + 



(14) 
(15) 



where -E para boiic an d {Api)\ mf , M correspond to the lower-order 
expansions |2} and |6|. 

The higher-order corrections have two major implications. 

Firstly, corrections to Api of quadratic or higher order in 
I lead to a deviation of the spectroscopic measurement [the 
right-hand side of Eq. |7])] from the actual average angular 
momentum (£). 

Secondly, corrections to the effective kinetic energy E of 
cubic or higher order in £ mean that even the very definition 
of the superfluid fraction, Eq. ([TJ, breaks down. A basic as- 
sumption of this definition is that the rotational properties of a 
(normal) gas are entirely characterised by its moment of iner- 
tia. This is correct provided the kinetic energy is a quadratic 
function of the angular momentum. Then, under a transforma- 
tion to a frame rotating at angular frequency u>, the interaction 
energy is unchanged, and the (parabolic) kinetic energy trans- 
forms as E -t E' = E + u) ■ (L) + \l d Lo 2 . Thus, for a 
state of given average angular momentum, the only material 
property characterising the net energy change is the moment 
of inertia 7 c i. For non-parabolic kinetic energy, however, the 
kinetic energy does not transform in any simple way, but de- 
pends on the populations of the individual angular momentum 
states and hence also on how interactions and/or temperature 
have distributed particles among these levels. The moment of 
inertia of even a normal gas cannot be assumed simply to be a 
constant I c \. 

We now turn to estimate the quantitative effects of these 
two forms of correction. 

First we consider the effects of the higher-order corrections 
to Api, Eq. (15) . Because of the departure of Apt from its 
linear expansion, we find Ap'£* ^ Apg- — Ap . This dif- 
ference is illustrated for different combinations of Vl R and At 
in Fig. |4] where the dashed lines show Ap'£* and the circles 
show the actual difference Ap^* — Apo. As a consequence, 
Eq. ([8J would give a systematically incorrect result even in the 
T — > limit. In practice this can be corrected by replacing the 
denominator in ([H) by Apt* — Apo. 

A more serious problem arises from the fact that non-zero 
temperature and/or interactions populate a range of £ states. 
Due to this broadening of the distribution function, the higher- 
order corrections to Ape in Eq. ( fT3] l lead to the situation that 
for a normal gas (or a relaxed superfluid) with average mo- 
mentum (£) = I* , the spectroscopic signal Ap is not just de- 
termined by the signal of the state at £*, but depends on the 
overall distribution function, so Ap ^ Apg- . Similarly, for 
the (metastable) superfluid with (£) = 0, one has Ap ^ Apo. 
Hence one would incorrectly determine the superfluid frac- 
tion. In order to estimate the size of this error, we here 



consider the effect of temperature using the example of an 
ideal Boltzmann gas. (Interactions will be considered in Sec- 
tions [V] and |VI]) We populate the different £ states in the low- 
est band (see Fig. [3]) according to a Boltzmann distribution 
n,£ oc exp[— ejj(£)/fcT]. The results are shown for a range of 
temperatures by the solid lines in Fig. |4] If the effects of ther- 
mal broadening were negligible, all these curves would agree 
with Api* — Apo (the circles in Fig. |4ji. 

Increasing il R decreases the relevance of higher-order cor- 
rections, but also diminishes the signal size Ap'i* and thus 
leads to a bigger experimental uncertainty. Increasing A£ has 
the opposite effect. One goal of this paper is to find parame- 
ters for which a good compromise can be reached, such that 
there is both a large experimental signal and small systematic 
inaccuracies. 

For 5 /fin "C 1, the influence of fl R and A£ on signal 
strength and higher-order corrections can be seen in the ex- 
pressions from perturbation theory. Using the perturbative ex- 
pansions in l/Cl R (cf. Section III i, the signal size is given by 



Ap'. 



Cl R MR 2 n R 



(A£) 2 + 0(l/n%), (16) 



and considering the corrections to Ap^ in (TT3J, the second- 
order coefficient d is given by 



-3^2 



8 ( h 



n R \MR 2 n R 



(A£) 2 + 0(l/n 4 R ). (17) 



These expressions show that increasing A£ is an effective way 
to increase the signal size, but also has an adverse effect on 
the accuracy. Accuracy is improved by increasing fl R , but at 
a reduction in signal size. As a trade-off one would thus try 
to make At as high as possible, and then increase fl R as long 
as the signal size remains big enough. As a reasonable, and 
experimentally feasible, compromise we choose At = 50 and 
fl R = 2ir x 100 kHz [as in Fig.|4](c)]. These are the values 
we will use through the remainder of this paper. 

Second, and finally, we return to the definition of super- 
fluidity in Eq. ([T}, and estimate its adequacy. Measuring the 
relative deviation ((t) —£*)/£* is a way to assess the depar- 
ture from parabolicity due to the terms of E which are odd 
in £. For the Boltzmann gas of Fig. |4] the relative deviation is 
maximal at £* ~ 6 — > (with a maximum value oc 1 and 
decreases with higher 6. At T = 1000 Hz [24|, the maximum 
value of the deviation at 6 — > is about 0.3%. To further 
quantify the size of non-parabolic corrections to the disper- 
sion relation, we have carried out fourth-order perturbation 
theory [25 1 in l/fl R . The next terms in the series in £ are 



E ~ -Eparabolic — 



h i (A£Y : 



1 ft 5 (Al) 4 4 
2^2 {MR 2 )±n 3 R 



0(£ 5 ), 



(18) 



where both coefficients have higher-order contributions of 
0(1/VL B ). The contributions of cubic and quartic term rel- 
ative to the parabolic term (h 2 t 2 /2M* R 2 ) are of order 10~ A £ 
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(b) Cl„ = 2n x 50 kHz, At = 10 




(c) n„ = 2vr x 100 kHz, At=50 
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FIG. 4: (Color online) Results for a normal gas: Ideal signal size 
Ap'l* and temperature-dependent deviation from this as a function 
of the detuning, 5/Qr. Data are for a classical Boltzmann gas at 
several values of Qr and At. Circles denote the zero-temperature 
limit (Ape* — Apo) at different values of £* . The curves are labelled 
from top to bottom, and £* is labelled from left to right. 



and 10 7 £ 2 , respectively. This is on the same order as the 
0.3% estimated above: For T = 1000 Hz and I* = 1, the root 
mean square angular momentum is yijP) ~ 16, so we would 
expect the relative correction due to the cubic term to be about 




FIG. 5: The system geometry employed in this section, a toroidal 
trap with hard walls and a cross-sectional area L x L. For the radius 
R 2> L, angular momentum £ can be unrolled unto the p z direction, 
p z = h£/R, which then has effective periodic boundary conditions. 
Thus £ and p z will be used interchangeably. 



0.2%. This deviation of the dispersion relation from parabol- 
icity sets an upper bound to the accuracy with which one can 
determine the superfluid fraction. 



V. INTERACTIONS 

In addition to the effects of non-zero temperature on a nor- 
mal fluid, we want to know what happens to Ap due to in- 
teractions in a superfluid. It is important to note that a gas is 
only superfluid because of interactions; in the limit of vanish- 
ing interaction strength, the critical velocity (see below, Sec- 
tion f 



V C i goes to zero, so a metastable superfluid flow cannot 



be maintained. 

Here we do not attempt a full model of the interacting gas 
in a trap. Rather, we aim to use a method that includes interac- 
tions in the simplest way. We therefore consider an interacting 
gas with uniform density. The density inhomogeneity near the 
walls is on the scale of the healing length £ = (Snna)" 1 ^ 2 , 
which is the distance over which the condensate recovers its 
bulk value from zero density at the walls. We assume that 
£ -C Ai? and that we can thus neglect the density inhomo- 
geneity. We therefore model the gas in a ring-like trap by 
considering a torus of radius R and with a rectangular cross- 
sectional area A — L x L, cf. Fig. [5] Since we neglect density 
inhomogeneities, we consider the system to be uniform and 
impose periodic boundary conditions in the transverse direc- 
tions. For L ~ Ai? <C R we can unroll the angular momen- 
tum from the circumferential direction of the torus onto a line 
with periodic boundary conditions, equivalent to a linear mo- 
mentum p z — (h/R)£, where £ is the angular momentum (in 
units of h) defined above. With this understanding we will use 
£ and p z interchangeably. 

We are only interested in the azimuthal motion, p z or £, in 
which we would like to measure superfluidity. Therefore we 
will integrate out the two transverse directions p±. These have 
a simple kinetic energy dispersion, p\ /2M. In the azimuthal 
direction we allow for a general dispersion relation ejj (p z ) , set 
by the energy of the lowest band of the Hamiltonian The 
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total dispersion relation is given by 



e°(p) 



1 

2M 



Pi 



(19) 



Later on we will compare the superfluid fraction with the con- 
densate fraction, so it is instructive to determine the number of 
excited particles 7Y ex . Here, we will determine this for an ideal 
non-interacting gas, as many of the features will carry over to 
the interacting case discussed below. In a non-interacting gas 
the number of excited particles is given by 



p#0 



(20) 



where ns(e) = [exp(e/fcT) — is the Bose-Einstein dis- 
tribution in the condensed system (for which the chemical po- 
tential vanishes, /i = 0). We split the sum, 



N ex = N e i( Pz =0)+ ]T N e i( Pz ), 



(21) 



defining 



n B 



(pi/2M + e|f( Pz )) . (22) 



Here we assume without loss of generality that the energy 
minimum is at ejj (p z = 0) = 0. We switch to an integral rep- 
resentation according to ^2 p — > J dp. The term N^(p z ) 
with p z 7^ does not pose any problems. However, when 
evaluating the p z = contribution, we find 



A 



J d 2 p ± n B (p 2 j2M) 

dppn B (p 2 /2M), (23) 



(2irh) 2 
2-kA 
~ (2^) 2 h 

which is infrared divergent. For p 0, the integrand is 
approximately given by pj (p 2 /2M kT) ~ 1/p, thus the in- 
definite integral at small momenta is ~ — log(p m j n ). Fortu- 
nately, in the thermodynamic limit this logarithmic divergence 
is unproblematic. Noting that the low-momentum cut-off is 
Pmin ~ h one finds N e i( Pz = 0) ~ (A/ X 2 T )\og(L/ X T ), 
where At is the thermal de Broglie wavelength. The other 
contributions are extensive, N^(p z ^ 0) ~ (V/X^), so 



N£(p z =0) At 

^ei(P^O) ~ R 



log(L/A T ), 



(24) 



and N^(p z — 0) can be neglected in the limit At <C R- The 
same reasoning applies when taking interactions into account, 
so from now on we will drop the p z = contribution to N ex 
and replace J2 P ^o b y E P ,#> Ep^o- 



A. Weak interactions with asymmetric dispersion 

In the following we consider an interacting gas, with a gen- 
eral, possibly asymmetric, non-interacting dispersion relation 



e°(p). Since we study atoms in the ultracold limit, all inter- 
actions can be treated as contact interactions, U e ff(r,r') = 
Uq5(t — r'). Here we take Uq = iith 2 a/AI, assuming for 
simplicity that the scattering length a is independent of the 
internal state. The Hamiltonian in the momentum basis is 

ijkl p.p' .q 

The indices k, I stand for internal states, e.g. one of the 
three hyperfine bands. As we assume that only the lowest 



band is occupied (cf. Section Ml, we can ignore inter-band 
mixing and simplify Eq. ( f25| to the Hamiltonian we will con- 
sider hereafter: 



H 



U 



z ^e-{p)a+a p + ^ d+ +q a+_ q d p >a p . (26) 
p,p',q 



I. Bogoliubov transformation 

We assume that the condensate has only one macroscop- 
ically occupied state p c = 0. Should the condensate be in 
a state p c ^ 0, we can shift all momenta by relabelling the 



states to p' = p — p c with energy e (p') 



e°(Pc 



P' 



We note that, for a superfluid, the state p c is not necessarily 
the ground state: for example, as in the protocol described 
above, a superfluid may remain condensed in £ = even af- 
ter a gauge field is applied such that the lowest-energy single 
particle state has £* ^ 0. Our discussion in this section also 
applies to these metastable superfluid states. 

The creation and annihilation operators of the ground state 
are d$ and So, respectively. We have (a^ do) = No, where 
iVo is the number of particles in the condensate. For No 
1, we can therefore neglect the commutator [&o,ad]~ = 1 
compared to the operator d^do, and replace the ground state 
operators by a c-number: clq ~ ~ V^Vo- Expanding the 
interaction term and only keeping terms of order O(Nq) or 
higher (i.e. assuming the only relevant interactions are with 
the condensate state), one finds 



H 



Ne°(0) 



UqN 2 
2V 



+ Z [e (p)-e°(Q) + U n ] a+a p 

p^O 

+ -g- 2^ {% a - P + a P a- p ) 

p^O 

i 

J2{[e°(p)-e°(0)+e 1 }a+a p 

p#0 

+ [e a (-p)~e°(0)+e 1 }d± p d_ p } 
/ 

+ ^ ei (a+at p + a p a_ p ) + E offset , (27) 

p#0 
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where the prime on the sum denotes that we only sum over 
(an arbitrary) half of momentum space. Here no = Nq/V is 
the condensate density, and we have defined ei = UqUq = 
AttH 2 ano / M , which sets the chemical potential, fi = e\. The 
constant terms have been absorbed into -Eoffset- 

We diagonalize the Hamiltonian using the Bogoliubov 
transformation 



with [d,d + ]_ 



— p — ffip ■ HtClf *p 

= [ / 8, J 8+]_ = 



-v 



(28) 



man 

,,2 



r = 1. The Hamilto- 

( p7| ) is diagonalized by choosing u 2 = h (- + l) and 
| (| - l), leading to 

1 



H = E, 



offset 



^ 2 

p#0 



E 

p#0 



(e + 7) a+a p , (29) 



where e eX c(p) = + 7(p) is me energy of Bogoliubov 
excitations. We have defined l26l 

1 

2 
1 

2 

e(p) = 7(p) + ei, 



7(p) 

70) 



[«°(p) 



} (-p)]-e°(0), 



e(p) = Je 2 - ef = v^7 + 2e0- 



(30) 



7(p) denotes the symmetric part of the non-interacting dis- 
persion relation, shifted such that j(p c ) — 0. e(p) is the sym- 
metric part of the Bogoliubov excitation energy. 7(p) denotes 
the asymmetric part of the excitation energy, and is the same 
for both non-interacting and interacting cases. As opposed to 
the standard case of a symmetric dispersion relation, we find 
7/O and a dependence of the energy e exc of Bogoliubov 
excitations on £* . We shall use this in Section [V CI to deter- 
mine the critical velocity at which the metastable superfiuid 
becomes unstable. 



2. Depletion of the condensate 
We can write the operator for the total number of particles 



as 

p#0 

When we apply the Bogoliubov transformation, we find 

1 re 



(31) 



iV^ + £ 2 (--l)+£-<d+d p ), (32) 

p#0 ' p^O 



where we have already taken the expectation value. At zero 
temperature, 



N(T = 0)=N Q + J2l(l-l) 



(33) 



After changing the sum to an integral, we can analytically in- 
tegrate out the transverse directions. In terms of the number 
density n = N/V we get 



n(T = 0)=n + 



n o a \ ■* 



1 



1_ 
(1 



l-*/l 



(34) 

where 7 is evaluated at p± = 0; 7 = 7(p 2 ;pj_ = 0). Each 
term of the sum gives the corresponding Estate occupation 
n e (T = 0). 

To find the temperature-dependent depletion, we calculate 
the expected excitation number: (dpd p ) = ub (e e xc(p))- The 
energy of such thermal excitations is given by e exc = e + 
7. Changing the integral over transverse directions into an 
integration over the dimensionless x = p 2 L /(2Me 1 ) leads to 



n(T) - n(T = 0) 



no a 



E 



dx- 



x + 1 + 7/ei 



1 



exp 



V / (a: + l + 7/ei) 2 -l + 7 /kT 



- 1 



(35) 



again with 7 = j(p z ',P± = 0). In general, this integral has 
to be evaluated numerically. Similarly to above, the terms 
of this sum give the finite-temperature distribution function 

n e (T)-n e (T = 0). 



B. Popov approximation 

To extend the validity of our approximation to higher tem- 
peratures, we employ the Popov approximation. In the homo- 
geneous system, the Popov approximation leads to the same 
Hamiltonian as in the Bogoliubov theory, except for a sim- 
ple excitation-independent energy offset. Thus the distribu- 
tion function stays the same — the only difference is that now 
we need to determine no self-consistently, both at zero and at 
non-zero temperature Il27ll28ll . 

As the zero-temperature distribution ([34} is symmetric 
around I = 0, at T = the average angular momentum is 
always zero, (I) = 0, for any interaction strength. This re- 
covers the expectation that p s /p = 1 even if the condensate 
is depleted by interactions. We calculate the zero-temperature 
depletion for a parabolic dispersion relation ejj = ^mb? = 7- 
We approximate the sum in Eq. ( |34| ) by an integral, 



n — no 
n 



a 2-kR 
2vri? 2irh 



/ dp z 




Jo 





1 



1 



(36) 



which can be evaluated analytically, giving 



n — no 
n 




p^O 



3^ 



(37) 
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where we substituted ej = 4Tth 2 ano/M. The resulting ex- 
pression for the condensate fraction, 



no 
n 



1 



,1/2 



(38) 



is not in a closed form; we have to find a self-consistent so- 
lution for n numerically. At small depletions n ex /n <C 1, 
Eq. ([38]l simplifies to the Bogoliubov result ll27l . no/n — 



1_ « („ ffl 3)l/2. 

According to this theory, one would simply occupy the low- 
est band as stated in the interacting distribution function given 
by Eq. ( |35j ) until convergence is achieved. For a parabolic 
dispersion relation, Eq. ( [35] > has an asymptotic behaviour 
n,£ oc l/£ 2 at large momenta, I ^> 1, and thus extends to very 
high £ However, the physical distribution of atoms does not 
extend beyond £ m . dx ~ R/R e , where i? e is the range of the in- 
teraction potential (in practical terms, this is of the same order 
as the scattering length a). Thus we do not trust the expression 
for m, Eq. ( |35] l, for l/R > 1 / R e ~ 1 /a, and so it is pointless 
to do the whole sum. In any case, all the important physics 
characterising the superfluid response is contained in p < h/£ 
(equivalent to £ < R/0- Hence we introduce a cut-off by ig- 
noring states with £ > £ cut . We choose £ cut ~ 10 x 2ttR/£ i , 
for which the relative error in Ap KSF — Ap SF due to the cut- 
off with respect to £ cut = oo is on the order of 10~ 4 . 



C. Critical velocity 

In order that the system is (meta)stable, the minimum 
of the excitation energy e exc = e + 7 has to be positive: 
ming e exc (£) > 0. This leads to a critical angular momen- 
tum shift £* A above which, for \£*\ > £* rit , the superfluid 
flow is unstable. (This is the Landau criterion for stable su- 
perfluid flow.) The determination of the critical £* in the gen- 
eral case requires a full (numerical) determination of the dis- 
persion relation of the lowest dressed band, ejj(p 2 ). In the 
limit in which the non-interacting dispersion relation can be 
taken to be parabolic [cf. Eq. Q], the symmetric and anti- 



symmetric parts of efl are given by -f(£) 

7W : 



A^£ 2 and 



h 2 r 



citation at given I is 



2M*fi 2< 

Then the lowest-energy Bogoliubov ex- 



£cxc(Rl = 0) = e + 7 



' 7i 2 ei 
M*R? 



% 2 t ^ 
M*R 2 ' 



(39) 



This is positive for all £ only if 

fi 2 \£*\/M*R 2 < yJh 2 e 1 /M*R 2 . 
Substituting e\ = AitTi 2 ano / M , this corresponds to 



m<Cit 



4irR 2 an 



M* 
W 



(40) 



If \£* I < £* clit , then e exc (£) > for all £. On the other hand, if 



and hence the system is unstable. It will relax due to the intro- 
duction of vortices. Eq. ( |40| shows that the metastable super- 
fluid will become unstable when the condensate density is too 
low (i.e. when the temperature is too high), when interactions 
are too weak, or when £* is too large. 



VI. RESULTS 

In order to present the results of our calculations for an 
interacting gas, we consider two scenarios. In the first one, 
we study a Bose-Einstein condensate (BEC) with weak inter- 



actions, Section VI A for which we show that the proposed 
measurement technique gives the expected feature that the su- 
perfluid density and condensate density are almost the same. 
In the second scenario, Section [VI B| we turn to the more in- 
teresting case of a BEC with strong interactions, for which the 
condensate fraction and superfluid fraction are significantly 
different from each other. We show that the measurement pro- 
tocol can distinguish the superfluid fraction from the conden- 
sate fraction. Finally, we will analyse the trade-off between 
signal strength and accuracy in Section fVIC| 

As described above, in our studies the condensate is as- 
sumed to be uniform and boundary effects are neglected. 
When discussing temperature dependences, we shall also as- 
sume that the density n is kept constant. For a harmonically 
trapped gas, with falling temperature the peak density quickly 
rises by more than an order of magnitude from near the transi- 
tion point to when most of the atoms are in the condensate. We 
ignore this effect, since we are primarily interested in the low 
temperature properties of the BEC and less in the behaviour 
around the transition point. In the results we present, we 
choose the following parameters: M = 23 m p , R = 10 pm, 
il R = 2tt x 100 kHz, A£ = 50, and n = 10 14 cm- 3 . These 
are typical parameters achievable with current technology. 



|^* I > £* r i t , then there always is some £ such that e c 



< 



A. Weakly interacting BEC 

For a new experimental technique, it is important to estab- 
lish that it works in a situation where the result is known. To 
that end, it will be valuable to see the results for a weakly 
interacting BEC, which can be well described by mean-field 
theories. Here we present the signal that one can expect to 
measure in such an experiment, for a weakly interacting BEC 
with nil 3 = 10~ 4 . The condensate fraction at T = is ap- 
proximately 99%, showing negligible depletion. 

As described in Section [II] the appearance of superfluidity 
displays a clear qualitative signal in this measurement tech- 
nique. The processes of cooling and of imposing non-zero 
gauge field £* do not "commute" for a superfluid, so its re- 
sponse is hysteretic. First introducing £* (by increasing the 
Raman detuning S) and then cooling the gas leads to a con- 
densate in the new ground state, £ c = £*. This is the "relaxed 
superfluid" (RSF), with spectroscopic signal Ap RSF . On the 
other hand, cooling to below the transition point at £* =5 = 
creates the superfluid in £ c = 0. Subsequently increasing 5 
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leads to a metastable (non-relaxed) superfluid (SF), with spec- 
troscopic signal Ap SF . Hence it is possible to measure two 
separate curves, Ap RSF respectively Ap SF . In contrast, for a 
normal fluid, cooling and increasing S do commute. 

In Fig. [6] (a), we show the spectroscopic signatures of 
the metastable superfluid (Ap SF ) and the relaxed superfluid 
(Ap RSF ) at temperatures below the BEC transition tempera- 
ture T c . (The raw signal has been scaled by £* so that the 
numbers are of the same order.) One can see that the curves 
of the metastable superfluid (Ap SF ) terminate at certain tem- 
peratures. This termination happens at a temperature (below 
the transition temperature) when the condensate becomes suf- 
ficiently depleted so that \£* | > £* liv Thus curves with higher 
\£*\ break off at lower temperatures. The difference in the 
signals Ap RSF ^ Ap SF below this temperature is a clear sig- 
nature of the metastable superfluid flow. 

In order to extract a quantitative measure of the superfluid 
fraction from these results, one could take the curves for the 
metastable superfluid Ap SF and use these in Eq. Appli- 
cation of Eq. §8§ requires knowledge of Apo and Ap' . As we 
have discussed, the result will also include some quantitative 
corrections from higher-order terms in the Taylor expansions 
of Eg and Api ( fl"5j ). In view of these facts, it is helpful to take 
the difference Ap RSF — Ap SF . This removes the dependence 
on Apo and also removes additive systematic errors. The dif- 
ferences Aj? RSF — Ap SF are shown in Fig.[6](b). Finally, we 
scale these differences by the value of the same quantity that 
is obtained at T — for a weakly interacting gas, 

Ap RSF (T = 0, na 3 « 0) - Ap SF (T = 0, na 3 « 0) . 

[Since here we consider na 3 — 10~ 4 , this is almost identi- 
cal to scaling by the T = limit of the curves in Fig. |6](b).] 
These measurements are for states in which we know the sys- 
tem should be completely condensed and perfectly superfluid 
at T = 0. They provide a direct measurement of the quantity 
Ap'£* , while also removing some multiplicative systematic 
errors. 

The final scaled curves are shown in Fig.|6](c). These show 
what the spectroscopic measurement would give for the super- 
fluid fraction of the weakly interacting BEC as a function of 
temperature. The different values of £* correspond to different 
effective rotation rates. For comparison, the dashed line shows 
the numerically determined condensate fraction no/n ||29ll . In 
this case, of a weakly interacting BEC, it is expected that the 
condensate and superfluid fractions should coincide. Fig.[6](c) 
shows that this result is recovered to a very good accuracy 
by the spectroscopic measurement of the superfluid fraction. 
Furthermore, in a practical experimental measurement with a 
harmonically trapped gas we expect to find even better agree- 
ment between the spectroscopic measurement and p s /p- here 
we chose n = 10 14 cm~ 3 , which is a representative value for 
the density in the BEC, but this leads to a transition temper- 
ature which is several times higher than that found in exper- 
iments. In a typical experiment, the density would vary such 
that T c is lower, so kT/h£lji would be smaller and quantita- 
tive corrections should have less of an effect. 
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FIG. 6: (Color online) Qualitative behaviour of the proposed su- 
perfluidity measurement as a function of temperature at na? = 
10" 4 : (a) the raw Ap signal scaled by (b) difference signal 
Ap RSF — Ap SF ; (c) difference signal normalised to 1 at T = 0, and 
comparison with the condensate fraction no/n. Curves break off due 
to reaching \£* | > £* A ~ V '4TrR 2 ano. The right edge of the graphs 
corresponds to the ideal gas transition temperature T c w 31.1 kHz. 
The curves are labelled from top to bottom. 
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B. Strongly interacting BEC 

We will now consider how the spectroscopic method can be 
used to provide an accurate quantitative measurement of the 
superfluid fraction in a regime where the condensate fraction 
no/n and the superfluid fraction are significantly different. 
To this end, we shall consider a relatively strongly interact- 
ing BEC, with an interaction parameter up to na 3 = 0.1. We 
cannot trust Popov theory to be an accurate quantitative theory 
of an atomic gas in regimes of strong interactions. However, 
we can still use the results of Popov theory to assess the accu- 
racy of the spectroscopic method of measuring the superfluid 
fraction. Specifically, we know that, even for a strongly inter- 
acting gas, in the low temperature limit the superfluid fraction 
should be p s /p = 1, while the condensate fraction can be 
significantly depleted. Here we wish to demonstrate that the 
difference between condensate fraction and superfluid frac- 
tion can be observed using the spectroscopic measurement 
method. 

As described in Section [rV] while considering an accurate 
quantitative measurement of p s /p, we face problems com- 
ing from higher-order corrections, which even cause the usual 
definition of superfluidity to break down. For this reason, in 
order to proceed we follow the protocol that was outlined in 



Section VIA We propose to define a spectroscopically mea 
sured superfluid fraction by 

Ap RSF (T, na 3 ) 
Ap RSF (T = 0,na 3 »0) 



P. 



spec 



p 



Ap SF (T, na 3 ) 
Ap SF (T = 0,na 3 



.0) 
(41) 

The denominator gives the reference signal at zero tempera- 
ture and with negligible interactions, where the system is com- 
pletely condensed into the ground state and 100% superfluid. 
(The interactions cannot be exactly zero, since one requires 
1 1* | < £* A ~ \/ATrR 2 an for metastability of the superfluid 



flow.) In principle, all four Ap values in Eq. (41 1 can be mea- 
sured in one experimental system by using a Feshbach reso- 
nance [30 1 to tune a. All measurements need to be taken at 
the same values for Hr, A£ and I* . One should take the limit 
I* — s- 0, analogous to lu — >• in the definition Q. 

As described in Section |IV| all higher-order corrections 
vanish for Hr — » oo, so in this limit the definition ([T| coin- 
cides with the spectroscopic definition ( |4Tj ). Thus, the super- 
fluid fraction is obtained from the spectroscopic signal ( |4T| ) by 
taking the limit 



Ps 

P 



lim 



spec 
P 



(42) 



Evaluating Eq. ( pT) at different values of £Ir gives a way 
to assess the influence of non-parabolic corrections. Further- 
more, this offers the possibility of improving the measurement 
by extrapolating to CIr — > oo. 

To analyse the quantitative accuracy of this method, we 
at first keep T = 0, i.e. we consider the accuracy of the 
method for the case of a perfect superfluid where we know 
that (£) — (or £* in the relaxed SF). We wish to show that 
at an interaction strength where the BEC is reasonably de- 
pleted, e.g. n ex /n ~ 0.1, the 100% superfluidity is still re- 
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FIG. 7: (Color online) Quantitative analysis of interaction effects 
at zero temperature: (a) raw Ap signal; (b) spectroscopically deter- 
mined superfluid fraction p s s pcc / p, and comparison with the conden- 
sate fraction no/n. Note that the condensate depletion is signifi- 
cantly higher than the deviation of p s s pcc /p from 1. The curves are 
labelled from top to bottom. 



covered by the spectroscopic method. The numerical results 
are shown in Fig. [7] as a function of increasing dimensionless 
interaction strength na 3 while keeping the density n fixed. 
Fig. [7] (a) shows the raw Ap results for the relaxed (Ap RSF ) 
and metastable (Ap SF ) superfluids. Fig. |7](b) shows the re- 
sulting spectroscopically determined superfluid fraction, from 



Eq. (41 1, using na = 10 to represent na 3 sw 0. For com- 
parison, in Fig [7] (b) the dashed line shows the numerically 
determined condensate fraction no/n for the SF OTI . The 
departure of the spectroscopic measurement of the superfluid 
fraction from 1 [the solid lines in Fig.|7](b)] shows that the er- 
ror in determining the superfluid fraction by the spectroscopic 
method increases with interaction strength. However, this er- 
ror is much smaller than the excited fraction n ex /n. Therefore, 
the spectroscopic method allows one to distinguish clearly be- 
tween the condensate fraction and the superfluid fraction in a 
strongly interacting BEC. 
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FIG. 8: (Color online) Trade-off between signal size and accuracy 
at zero temperature. Signal size is assumed to be given by Ap'l*; 
compare Fig.|4](c). Accuracy is determined by the deviation of the 
spectroscopic measurement l |4 1 ^ from the expected unity. This devia- 
tion has been calculated at na? = 10 -2 (corresponding to a conden- 
sate depletion n cx /n m 12%), amounting to a vertical slice through 
Fig. 0(b) atna 3 = 10" 2 . 



FIG. 9: (Color online) Comparison of p s s pcc / p and p 3 /p = 1 — 
(I) /£* at na 3 — 10 -2 as a function of temperature. As in Fig. [fl] the 
curves break off for \£*\ > l* tit ~ V 'AirFPano. The right edge of 
the graph corresponds to the ideal gas transition temperature T c w 
31.1 kHz. The curves are labelled from top to bottom at the left-hand 
edge of the graph. Each curve starts at the y-value corresponding to 
the one in Fig.[7](b) at na 3 = 10" 2 . 



C. Optimal experimental parameters 

In the preceding sections we have shown that the spectro- 
scopic method is capable of providing both a qualitative sig- 
nature of superfluidity and a quantitatively accurate way to 
measure the superfluid fraction. We now turn to discuss the 
experimental parameters to use to optimise this technique as a 
quantitative measurement of superfluid fraction. 

As was discussed in Section|lV] there is a trade-off between 
the signal size and the quantitative accuracy of the method. 
In particular, we showed that it is advantageous to have as 
large a value of A£ as possible, and then to increase CIr as 
much as possible to improve accuracy but not so much as to 
make the signal Ap too small. From these considerations we 
were led to choose A£ — 50 (from practical considerations of 
achieving beams of high angular momentum), and we settled 
on a compromise value of Q# = 27r x 100 kHz. However, 
there still remains the question of what is the best detuning S 
(and therefore value of £*) to use. Formally, following Eq. <|8j 
we should consider the limit £* — > 0, such that the superfluid 
remains metastable even very close to T c . But in this limit the 
signal becomes very small. What is a reasonable value of £* 
to use? 

To explore this question, in Fig. [8] we present the results 
of calculations at zero temperature and na 3, = 10~ 2 , which 
corresponds to a condensate depletion of about 12%. 

In Fig. [8] the diamonds show the signal strength, as given 
by Ap'£* . This is the ideal difference between the signals of 
a normal fluid and a superfluid at T = 0, neglecting higher- 
order corrections to Ape. As in Fig.|4](c), this is a good ap- 
proximation to the actual signal. These results show that, for 
the parameters chosen, the fractional imbalance Ap must be 
measured to an accuracy of about 0.01 to 0.05. While this 



is a small fractional imbalance, it should be feasible to detect 
signals of this size by averaging over many shots. 

The circles in Fig. |8]display the inaccuracy in the spectro- 
scopic measurement of the superfluid fraction, i.e. the rel- 
ative deviation of p s s pec /p from the expected value of 1 at 
T = 0. Hence the circles in Fig. [8] correspond to a vertical 
slice through Fig.|7](b) at na 3 — 10 -2 . Here, an inaccuracy 
of 0.03 means that a 100% superfluid would be misinterpreted 
as only having 97% superfluid fraction. In contrast, recall that 
the condensate fraction is about 88%. So this still allows a 
clear distinction between superfluid and condensate fraction. 

The results in Fig. [8] would suggest picking \£*\ s» 35, 
where the inaccuracy passes through zero. However, one 
should note that these results only show the inaccuracy at 
zero temperature. In Fig. [9] we plot the ratio p s s pcc /p s as a 
function of temperature. Here p s s pcc is the spectroscopically 
measured superfluid density and p s is the expected superfluid 
density, as computed from the usual definition of superflu- 
idity (flj, p s /p = 1 — (£)/£*. The ratio stays close to one 
over a wide temperature range, showing that the spectroscopic 
method provides a good measure of the superfluid fraction. 

Combining the issue of signal size and accuracy over a 
range of temperatures, we find that (for the parameters stud- 
ied), a good choice is \£* \ between 20 and 30. In this range, 
the signal size is 0.04 — 0.05; the spectroscopically measured 
superfluid fraction is accurate to about 0.02 (much less than 
the depleted fraction of 0.12); the accuracy remains at this 
level up to ~ 0.9T C . 

VII. SUMMARY 

In summary, we have assessed the feasibility of a recent 
proposal ifTUl to measure the superfluid fraction of atomic 
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BECs. This proposal involves the use of optically induced 
gauge fields to simulate rotation, and allows the superfluid re- 
sponse to appear in a spectroscopic signal. We have calculated 
the expected spectroscopic signatures for three-dimensional 
BECs with uniform density. One conclusion of our studies 
is the demonstration that this technique can be used to ob- 
tain a qualitative experimental signature of the superfluid re- 
sponse (by comparing the spectroscopic signals when the or- 
der of turning on the gauge field and cooling the gas is re- 
versed), with a spectroscopic signal that is large enough to be 
detectable in experiment. Furthermore, and most importantly, 
we have shown that the technique can be used to extract a 
quantitative measurement of the superfluid density. The ac- 
curacy of the measurement technique can be improved at the 
expense of the size of the signal. Using realistic values for the 
experimental parameters, we have shown that a compromise 
can be reached where both (i) the signal is sufficiently large 



to be feasible to measure and (ii) the superfluid density is de- 
termined to sufficient accuracy to allow quantitatively useful 
information to be extracted. Notably, our results show that the 
technique can allow a clear experimental measurement of the 
distinction between condensate and superfluid fraction of a 
strongly interacting Bose gas. Our results support the useful- 
ness of this technique 1 10 1 for measuring the superfluid frac- 
tion, and provide guidance for the parameters required in fu- 
ture experimental implementations. 



Acknowledgments 

We are grateful to Jean Dalibard for helpful comments on 
the manuscript. This work was supported by EPSRC Grant 
No. EP/I010580/1. 



[1] M. H. Anderson etal, Science 269, 198 (1995). 

[2] C. Raman et al, Phys. Rev. Lett. 83, 2502 (1999). 

[3] K. W. Madison, F. Chevy, W. Wohlleben, and J. Dalibard, Phys. 

Rev. Lett. 84, 806 (2000). 
[4] C. Ryu et al, Phys. Rev. Lett. 99, 260401 (2007). 
[5] I. Bloch, J. Dalibard, and W. Zwerger, Rev. Mod. Phys. 80, 885 

(2008). 

[6] A. J. Leggett, Rev. Mod. Phys. 71, S318 (1999). 
[7] E. Andronikashvili, J. Phys. USSR 10, 201 (194 6). 
[8] S. Riedl et al, ArXiv e-prints, [arXTv:0907.3814| /2 (2009). 
[9] T. Ho and Q. Zhou, Nature Physics 6, 131 (2009). 
[10] N. R. Cooper and Z. Hadzibabic, Phys. Rev. Lett. 104, 030401 
(2010). 

[11] J. Dalibard, F Gerbier, G. Juzeliunas, and P. Ohberg, ArXiv 

e-prints. liXiv: 1008.5378V 1 (2010). 
[12] Bose-Einstein Condensation, edited by A. Griffin, D. W. Snoke, 

and S. Stringari (Cambridge University Press, Cambridge, 

1996). 

[13] L. Tisza, Nature 141, 913 (1938). 

[14] L. Tisza, Phys. Rev. 72, 838 (1947). 

[15] L. D. Landau, J. Phys. USSR 5, 71 (1941). 

[16] S. Franke-Arnold, L. Allen, and M. Padgett, Laser & Photonics 

Reviews 2, 299 (2008). 
[17] J. M. Blatt and S. T. Butler, Phys. Rev. 100, 476 (1955). 
[18] F. Zambelli and S. Stringari, Phys. Rev. Lett. 81, 1754 (1998). 
[19] F. Chevy, K. W. Madison, and J. Dalibard, Phys. Rev. Lett. 85, 

2223 (2000). 

[20] G. Juzeliunas, P. Ohberg, J. Ruseckas, and A. Klein, Phys. Rev. 
A 71, 053614 (2005). 



[21] I. B. Spielman, Phys. Rev. A 79, 063613 (2009). 

[22] Y.-J. Lin etal, Phys. Rev. Lett. 102, 130401 (2009). 

[23] For simplicity, we neglect any quadratic Zeeman effect, assum- 
ing that the magnetic field is small. 

[24] For easy comparison with other parameters, for example the 
Rabi frequency, we express temperatures in frequency units, 
such that 1 kHz = 48 nK. 

[25] F. M. Fernandez, Introduction to perturbation theory in quan- 
tum mechanics (CRC Press, Boca Raton, FL, 2001). 

[26] Since we assumed that the condensate is in p c = 0, the last 
term of "/(p) in i30i refers to the condensate state; in general 
the argument to e u has to be changed appropriately: "f(p) = 
|[f (Pc+p) + e (p c -p)]- e °(p c ). 

[27] C. J. Pethick and H. Smith, Bose-Einstein Condensation in Di- 
lute Gases; 2nd ed. (Cambridge University Press, Cambridge, 
2008). 

[28] V. N. Popov, Functional Integrals and Collective Excitations 
(Cambridge University Press, Cambridge, 1988). 

[29] The condensate fraction no /n deviates from the 1 — (T/T c ) 3 / 2 
behaviour expected in an ideal Bose gas due to weak interac- 
tions as considered within the Popov approximation. 

[30] C. Chin, R. Grimm, P. Julienne, and E. Tiesinga, Rev. Mod. 
Phys. 82, 1225 (2010). 

[31] For the parameters used, this is indistinguishable from the 
condensate density determined by self-consistently solving 
Eq. < |3 8[ > for no, showing that non-parabolic corrections are ir- 
relevant when considering the condensate fraction. 



